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ABSTRACT 

Context. X-ray studies of AGN with powerful nuclear winds are important to constrain the physics of the inner accretion/ejection 
flow around SMBH, and to understand the impact of such winds on the AGN environment. 

Aims. Our main scientific goal is to constrain the properties of a variable outflowing absorber that is thought to be launched near the 
SMBH of the mini-BAL QSO PG 1 126-041 using a multi-epoch observational campaign performed with XMM-Newton. 

Methods. We performed temporally resolved X-ray spectroscopy and simultaneous UV and X-ray photometry on the most complete 
set of observations and on the deepest X-ray exposure of a mini-BAL QSO to date. 

Results. We found complex X-ray spectral variability on time scales of both months and hours, best reproduced by means of variable 
massive ionized absorbers along the line of sight. As a consequence, the observed optical-to-X-ray spectral index is found to be 
variable with time. In the highest signal-to-noise observation we detected highly ionized X-ray absorbing material outflowing much 
faster ( vx ~ 16 500 km s _1 ) than the UV absorbing one (u,„ ~ 5 000 km s _1 ). This highly ionized absorber is found to be variable on 
very short (a few kiloseconds) time scales. 

Conclusions. Our findings are qualitatively consistent with line driven accretion disk winds scenarios. Our observations have opened 
the time-resolved X-ray spectral analysis field for mini-BAL QSOs; only with future deep studies will we be able to map the dynamics 
of the inner flow and understand the physics of AGN winds and their impact on the environment. 

Key words. Accretion, accretion disks - Methods: data analysis - Techniques: spectroscopic, photometric - Galaxies: quasars: indi- 
vidual: PG 1126-041 - X-rays: individuals: PG 1126-041 


1. Introduction 

The structure of the inner regions (sub-parsec scale) of Active 
Galactic Nuclei (AGN), as probed by UV and X-ray observa- 
tions, seems to be very complex and certainly it is still not well 
understood. The optical/UV continuum emission is most proba- 
bly due to the thermal emissi on from an optically thic k, geomet- 
rically thin accretion disk (IShakura & Sunvaevlll973l) . The disk 
surrounds the central supermassive black hole (SMBH), span- 
ning radii from a few up to several hundreds of gravitational radii 
( r„ = GMbh/c 2 ), i.e. possibly from the innermost stable cir- 
cular orbit around the SMBH up to the disk self-fragmentation 
radius (e.g., lPringlelfl98 It) , The origin of the X-ray continuum 
emission is less clearly understood and is thought to be the 
result of the comptonization of accretion disk UV seed pho- 
tons into a “cloud” of very hot electrons, the so-called X-ray 
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corona dHaardt & Maraschilll 991t lHaardt et alJfT994li . What is 
clear both from variability and microlensing studies is that the 
X-ray emission region is much smaller than the UV one and 
spans only a few up to tens gravitational rad ii (e.g.,|Mprgan et al. 
120081 : lOiartas et al.ll2009bt lDai et al]l2010h . An extensive X-ray 
monitoring of the inner regions of the nearby Seyfert 1.8 NGC 
1365 made possible to constrain the size of the X-ray emitting 
region to be again a few r g , using a totally independent technique 
(IRisaliti et al. 120093) . 

The UV and X-ray continuum photons are observed to be 
reprocessed (absorbed, re-emitted, scattered, reflected) by gas 
and dust in the inner regions of AGN, so that we almost never 
observe just the primary continuum emiss ion of these objec ts 
(see e.g. the comprehensive analysis of IRisaliti et al.ll2005hllaL 
120071 l2009al again on NGC 1365). Conversely, observing the 
reprocessing features in the UV and X-ray bands can put strong 
contraints on the geometry, the physical characteristics, and the 
dynamics of the gas hosted in the inner regions of AGN. 

In the last few decades, spectroscopy in the UV and X-ray 
band have revealed the presence of substantial column densities 
of ionized gas outflowing from the inner regions of AGN. 
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In the UV band we observe blueshifted spectral absorption 
lines due to resonant transitions of ionized metals such as Mg II, 
Al HI, Si IV, C IV, N V, O VI. Depending on the width of 
the absorption troughs, quasars hosting such features are clas- 
sified as Broad Abs orption Line Quasars (BAL QSOs, FWHM 
> 2000 km s" 1 . e.g.lTurnshek et al.ll 1 98(j IWevmann et al .ll 1 99 1 fc 
iReichard et al.l 12003 ). mini-Broad Absorption Line Quasars 
(mini-BAL QSOs, 500 km s 1 < FWHM < 2000 km s '), 
and Narrow Absorption Line Quasars (NAL QSOs, FWHM 
< 50 0 km s _1 ; see IGangulv et al.l 1200 it lHamann & Sabral 
120041) . BALs featu res are observed in ~ 15% of optically se- 
lected QSOs (e.g., iHewett & Foltz [2003b IReichard et al.l [20031 


iKnigge eT al. 2008; Alle n et al.l201 ll) . Mini-BALs and NALs are 
together observed in ~ _J 2 - 30% of optically selected QSOs 
IGangulv & Brothertonl2008i) . 

Despite the very different absorption trough widths 
observed, BAL, mini-BAL, and NAL QSOs share the 
same velocity ranges, reaching UV terminal velocities v\ 


uv 


from 


a few 1 (r km s 1 u p to several 10 4 k m s 1 
Rodriguez Hidalgo e t al.l 120071; iMisawa et al.l 12007 at 


(e.g. . __ 

Gibs on et al.l 2009hl) . Furthermore. is known to corre- 
late with the continuum lum inosity L uv (lLaor & Brandt! [2002: 
IGangulv & Brothertonll2008l ). meaning that the highest terminal 
velocities are observed in the objects with the highest UV con- 
tinuum luminosity, suggesting an important role of the AGN UV 
radiation pressure in the acceleration of such winds. 

In the X-ray band absorption due to ionized species, such as 
N VI- VII, O VII- VIII, Mg XI-XII, Al XII-XIII, Si XIII-XVI, as 


well as L-shell transitions of Fe XVII-XXIV, are observed to be 
blueshifted by a few hundred to a few thousand km s_ l in ~ 50% 
of type 1 AGN (the “warm absorber”, e.g. lRevnoldslll997t 
iPiconcelli et a Ll 120051 iMcKernan et al.ll200% . Moreover, in re- 
cent years, thanks to the high collecting area of X-ray satellites 
such as XMM -Newton, Chandra, and Suzaku, blueshifted ab- 
sorption lines due to highly ionized gas (i.e. Fe XXV, Fe XXVI) 
outflowing at much higher velocity ( v out ~0. 05-0. 2c) have been 
obse rved in a number of AGN, at both low and high redshift (see 
e.g. lPounds et al .1120031 ICannil2006l ITbmbesi et al.ll2010allbL and 
references therein). In particular. X-ray BALs in He-like iron, 
blueshifted up to u out = 0.7c, have been observed in gravita- 
tional ly lense d BAL and mini-BAL QSOs (IChartas et al.ll2002L 
120031 12009di . 


Theoretically, powerful winds can be launched from AGN 
accretion disk either by thermal, magneti c or radiative pres- 
sure (see e.g. iKonigll I2006L lEverettl l2007t IProgal 120071 for a 
review). The thermal mechanism can account for low velocity 
outflows arising at large r adii (e.g. the X-ray warm absorber, 
IChelouche & Netzeill2005 ); instead, in order to explain the ex- 
treme velocities observed in the other absorbers, either magnetic 
or radiation pressure must be invoked. For instance, line driven 
accretion disk wind models show that at high accretion rates 
the UV radiation pressure from the disk can be high enough to 
launch powerful winds from distances of the order of a few tens 
of r g . In the inner part of this flow high column densities of X-ray 
absorbing gas shield the wind from the strong X-ray continuum 
source, preventing it from becoming overionized (iMurrav et al.l 
119951 [Proga et al.ll2000t iProga & K allman 2004). Thus, the ge- 
ometry and the dynamics of such winds depend critically on the 
UV/X-ray flux ratio and on the AGN accretion rate. On the other 
hand, for the launch and acceleration of magnetically driven disk 
winds, there is no need for X-ray shielding, and the accretion 
rate is not a critical parameter; the ionization state of the gas, 
that i s a critical parameter, will however still depend on them 
(e.g JKonigl & Karti ell 1 9941 lEverettll2005 ). 


BAL and mini-BAL QSOs are actually known to be “X-ray 
weak” with respect to the average QSOs. A comparison of the X- 
ray fluxes of BAL and non-BAL QSOs with the same UV flux 
indicates that BAL QSOs typically have X-ray fluxes that are 
10-30 times lower than those of non-BAL QSOs (iGreen et ah 


1995 


2000; 


LaoretaQ 19971). A number of studies (e.g. iBrandt etal 


Gree n et al.ll200ll lGallagher et al. 2006) strongly support 


absorption to be the reason for the X-ray weakness. Recently a 
number of “X-ray bright” BAL and mini-BAL QSOs has been 
discovered, with the latter clas s being X-ray brighter and less X- 
ray absorbed than the former (iGiustini et al.ll200'8l iGibson et al.l 
l2009al) . NAL QSOs have been shown to be generally unabsorbed 
in X-rays^ and X-ray brighter than both mini-BAL and BAL 
QSOs (IChartas et alj|2009al) . 

The physical link between the UV and X-ray outflowing ab- 
sorbers, and the dynamical behavior of the wind, are far from 
being understood. Reported estimates of the distances of the ab- 
sorbers from the central SMBH range from the inner regions of 
the accretio n disk (e.g ..ICharta s et al.l20 09c). to the parsec-scale 
torus (e.g., Blustin et all 20051) . to the kilonarsec scales (e.g., 
iDunn et al.ll2010l) . Currently the uncertainties of the launching 
radii are relatively large, which translates to a high uncertainty 
on the mass outflow rate and on the kinetic energy injection 
associated with such winds. It could also be that the launch- 
ing radii depend on the mass of the SMBH. Constraining the 
physical mechanism responsible for launching and accelerating 
AGN winds, i.e. understanding the link between the accretion 
and ejection processes in AGN, can help not only in understand- 
ing the evolution of these sources, but also in quantifying their 
impact on the s urrounding environment, i.e. the amount of feed- 
back (see e .g. iDi Matteo et al.l 120051: IScannanieco et al.l [2005; 
lElvislUOObh . 


In this article we present the results of an XMM -Newton 
observational campaign on the mini-BAL QSO PG 1126-041, 
that provides both the largest dataset (four pointings) and the 
deepest X-ray exposure (130 ks) on a mini-BAL QSO to date. 
The source is extremey interesting in showing all the three ab- 
sorbers discussed above: mini-BALs in the UV, an X-ray warm 
absorber, and a highly ionized, high velocity X-ray absorbing 
outflow. Constraining the X-ray/UV properties of this source can 
give precious insights into the physics of AGN winds. 

The structure of the article is as follows: Section [2] sum- 
marizes the source properties known from literature; Section [3] 
reports on the new XMM -Newton observations and data re- 
duction. EPIC pn timing analysis results are first presented in 
Section [4l EPIC pn and MOS spectral analysis is then presented 
in Section[5](average spectra of each epoch) and[6](time resolved 
2009 Long Lookpn spectra). Section[7]contains the results of si- 
multaneous optical/X-ray photometric analysis. Results are dis- 
cussed in Section[8]and conclusions are presented in Section[9] 

A cosmology with Ho = 70 km s-1 Mpc' 1 , Oa = 0.73 and 
Dm = 0.27 dKomatsu et al.lf20TTh is used throughout the paper. 


2. PG 1126-041 

PG 1 126-041 (a.k.a. Mrk 1298) is a nearby (z = 0.06) radio quiet 
AGN with a luminosity in between tho se typical of Seyferts and 
QSOs (M b = -22 . 8. G.„i ~ 10 I2 U . [S~chmidt & Gree^l 119831: 
ISanders et al.ll 1 980L ISurace et al.l l200l l). Its optical strong F e II 
and weak [O III] line emission (iBoroson & Greeiil 1 1 992l) are 
characteristics of Narrow Line Seyfert 1 galaxies (NLS1), al- 
though the FWHM of its H/j line of 2150 km s' 1 is slightly 
larger than the 2000 km s ~* value used in defining NLSls 
(iShuder & Osterbrocklll98ll) . Its black hole mass is estimated 
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Table 1 . Main properties of PG 1 126-041 


RA 

1 lh 29m 16.6s 

(a) 

DEC 

-04h 24m 08s 

(a) 

Z 

0.062 

(a) 

D l 

284 Mpc 


N gaI 

iV H 

4.35 x 10 20 cm -2 

(b) 

E(B - V) 

0.055 mag 

(c) 

Mbh 

5.0 - 12 x 10 7 M o 

(d-e) 

m B 

14.92 mag 

(a) 

Lbol 

8.9 x 10 U LO 

(f) 

L 2-10 

2 x 10 43 erg s _1 

(h) 

j/ uv 

^out 

5000 km s~* 

(g) 

v x 

u out 

16500 km s~‘ 

(h) 


References: (a )IS chmidt & Greenl dl 983ll: IbllKal berla et ali d2005h : (c) 


(d) Vesten giard &_Peterson|j 20061): (e ) 

Ifl lSurace et al.l d2001h : (g) lWang et alj d 1999b : 


S chlecel et al 


Dasyra et ajj d200' 


(h) This work, referred to the 2009 Long Look observation. 


to be log(M BH /M 0 ) ~ 7.7 using the L,j(5100A) continuum lu- 
minosity density and the FWHM(HyS) measured from a single- 
epoch optical spectrum, and the empirical relations reported 
by IVestergaard & Petersonl ( 2006 ). However, the stellar veloc- 
ity dispersion in the host galaxy of PG 1 126-041 measured from 
the stellar CO absorption line width in the VLT H-band spec- 


tra gives_a slightl 
8.1 (Dasyra et al 
“soft X-ray weak' 


higher black hole mass of log(MBH/AL 0 ) 

PG 1126-041 is also classified as 


IM 

•0 by B 


Brandt et al.l d2000l) . who measured an 
a ox ~ -2.19 based on ROSAT X -ray observations and o n 
the optical flux density as given by iNeugebau er et alJ d 1 9871) . 
Tananbau m et alJ (119861) measured an a ox — -1.90 referring 
to an earlier Einstein X-ray observation and the 250 0 A rest- 
frame flux density as given bv ISchmidt & Greenldl983l) . Neither 
of these a ox measurements were based on simultaneous X-ray 
and optical/UV observations. IWang et al.l (1 9991 ) report on the 
IUE and ROSAT observations of PG 1126-041. The 1992 IUE 
observation shows a high velocity UV absorbing outflow as 
mini-BALs in the C IV, N V, and Si IV ions blueshifted up to 
~ 5000 km s 1 . The simultaneous ROSAT pointing revealed a 
0.2-2 keV spectrum best fitted by a steep powerlaw emission 
with F ~ 2.8, absorbed by a column density Nw ~ 3 x 10 22 cm -2 
of ionized gas, log£ ~ 1.7 erg cm s _I . A subsequent IUE spec- 
trum taken in 1995 showed a contin uum brighter by_a factor of 
two and weaker absorption lines. iKomossa & Meerschweincheni 
J200d) confirmed the detection of ionized absorption in the 1992 
ROSAT spectrum and also suggested X-ray temporal variability 
on timescales as short as 800 s. 

The main properties of PG 1126-041 are summarized in 
Table □ 


3. Observations and Initial Data reduction 

XMM-Newton observed PG 1126-041 four times, once in 2004 
December, twice in 2008 June and December, and once more 
in 2009 June. In Table [2] we list the main parameters of each 
dataset. We reduced all the datasets using the XMM-Newton S AS 
v. 10.0.0 and calibration files generated in August 2010. For each 
dataset we extracted the light curve from the whole field of view 
of events with energies greater than 10 keV for each European 


1 Soft X-ray weak sources are those being characterized by an ob- 
served optical-to-X-ray spectral index a ox < -2, see Section[7] 


Photon Imaging Camera (EPIC) MOS cameras, and with ener- 
gies 10 keV < E < 12 keV for the EPIC-pn instrument. We 
then filtered out the time periods in which there was signifi- 
cant background flaring. Source counts from the cleaned event 
files were extracted from circular regions centered on the PG 
1 126-041 position, with ~ 3 0-45” radii, depending on the signal- 
to-noise ratio (S/N) as determined with the eregionanalyse 
task. Background counts were extracted from circular regions 
of the same area as the source ones. We retained up to double 
events (flagged as #XMMEA_EP) for the pn dataset and quadru- 
ple events (flagged as #XMMEA_EM ) for the MOS datasets. 
For each dataset we generated the ancillary response file and 
the redistribution matrix file at the source position using the 
tasks arfgen and rmfgen. Source light curves were created 
in different energy bands and corrected with the epiclccorr 
task for both the local background and the mirrors/detector in- 
efficiencies. Due to the low X-ray fluxes involved, pile-up ef- 
fects are negligible in our spectral analysis. For the same rea- 
son we cannot make use of the reflection grating spectrometer 
(RGS) aboard XMM-Newton. Optical Monitor (OM) data were 
all taken in Image mode and were processed using the standard 
procedure as suggested by the XMM-Newton SOC. 


4. Timing Analysis 

The left panel of Fig. Q] shows the background-subtracted 0.2- 
10 keV pn light curves of PG 1126-041, for each epoch of 
observation, binned to 1 ks. Compared to the light travel time 
ti = r g /c, the time bin corresponds to ~ 1.5 - 4 r g for the two 
different estimates for the black hole mass of PG 1 126-041. 

The variations on month timescales are dramatic: the source 
flux increases by a factor of ~ 4 from December 2004 to 
December 2008, then comes back to a low count-rate regime in 
June 2009. For comparison, the average 0.2-10 keV pn spectra of 
each epoch of observation are shown in the right panel of Fig. Q] 
Most of the spectral variability on month time scales occurs at 
energies E < 6 keV. The spectral shape is peculiar and obviously 
deviates from a simple power law. A prominent broad absorption 
feature is evident in all the spectra around E ~0.6- 1 .5 keV. 

The high signal to noise (S/N) ratio and the long contigu- 
ous good exposure time allowed us to perform a more detailed 
timing analysis on the 2009 Long Look pn exposure. 

Figure [2] summarizes the results. The left panel shows the 
Root Mean Square variability (RMS; Vaug han et al.l 120031: 
iPonti et akl 120041) computed on the 0.2-10 keV band: on very 
short time scales (ks) the source flux is varying up to 50% at 
energies E > 1.5 keV, while it is constant at lower energies. 
The middle panel shows the background subtracted lightcurves 
extracted in the 0.2-1. 5 and 1.5-10 keV band; following the 1.5- 
10 keV band flux fluctuations, we splitted the exposure in eight 
time intervals on which we perform temporally resolved spectral 
analysis. Three representative spectra extracted during the 2nd, 
the 7th, and the 8th time interval are shown in the right panel: 
most of the spectral variability on ks time scales occurs at 1 .5 
keV < E <9 keV. 

From the timing analysis, it follows that the X-ray spectra 
of PG 1126-041 can be divided in two spectral components: 
one dominates at £ > 1.5 keV and contributes 100% to the ks 
time scale variability, while the other one dominates at lower 
energies and is significantly variable only over longer (month) 
time scales. With these considerations in mind, we proceed with 
the analysis of both the average (Sect. [5j and the time resolved 
(Sect. [6J spectra. 
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Table 2. X-ray observation log 


OBSID 

Date 

Duration 

Exposure 

Mode 

Filter 

Count-rate 

l°g fo.2-2 

log h- to 

Epoch 




Ml /M2/pn 



MOS/pn 

pn 

pn 




[ks] 

[ks] 



[0.1 ct s" 1 ] 

[erg cm -2 s -1 ] 

[erg cm -2 s -1 ] 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

0202060201 

12/31/2004 

33.8 

32.9/32.9/28.3 

LW 

t 

0.61±0.02/ 2.01 ±0.03 

-12.89 

-11.93 

Dec. 2004 

0556230701 

06/15/2008 

31.4 

3.7/ 3.5/ 3.2 

FF 

m 

0.91±0.05/3.01±0.10 

-12.57 

-11.76 

Jun. 2008 

0556231201 

12/13/2008 

11.9 

4.9/ 4.7/ 3.7 

FF 

m 

1.68±0.06 / 5.93±0.20 

-12.31 

-11.60 

Dec. 2008 

0606150101 

06/21/2009 

134.3 

85.4/89.2/91.2 

FF 

m 

0.57±0.01 / 1.96±0.02 

-12.81 

-11.99 

2009 Long Look 


Notes. Col.(l): Observation ID; Col.(2): Date of observation; Col. (3): Nominal duration of the observation; Col. (4): Net exposure time for each 
instrument after the background flaring filtering was applied; Col. (5): Observing Mode, LW= Large Window. FF=Full Frame; Col.(6): Optical 
Filter applied, t=thin, m=medium; Col. (7): Net count-rate in the 0.2-10 keV and 0.3-8 keV range for the EPIC-pn and EPIC-MOS instruments, 
respectively, after the local background subtraction; Col. (8): Observed EPIC-pn flux in the 0.2-2 keV band; Col. (9): Observed EPIC-pn flux in the 
2-10 keV band; Col. (10): Epoch name as used in the text. 
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C 

3 

C 
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Fig. 1. Variability of PG 1 126-041 over month time scales. Left panel: background subtracted 0.2-10 keV pn light curves extracted in 
the four different epochs of XMM-Newton observations and binned to 1 ks: December 2004 (blue squares), June 2008 (red circles), 
December 2008 (green triangles), and 2009 Long Look (black circles). Right panel: the time-averaged pn spectra corresponding to 
the time intervals shown on the left and using the same color codes. 
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Fig. 2. Variability over kilosecond time scales detected during the 2009 Long Look observation. Left panel: RMS variability com- 
puted for the pn observation over the 0.2-10 keV band; arrows mark upper limits, error bars are at lcr confidence level. Middle 
panel: 0.2- 1.5 and 1.5-10 keV background subtracted pn light curves binned to 2 ks; the numbers in the 1.5-10 keV light curve mark 
the eight time intervals on which time resolved spectral analysis was performed. Right Panel: three representative spectra extracted 
during the 2nd, the 7th, and the 8th time intervals of the 2009 Long Look observation. 


5. Spectral Analysis: the average spectra 

In our spectral analysis we used the Xspec v. 12.6.0 soft- 
ware. All the spectra were grouped to a minimum number of 
20 counts per energy bin in order to apply the^ 2 statistics in the 
search for the best fit model (lAvnilll9'76h . Errors are quoted at 
lcr confidence, unless otherwise stated . All the models include 


Gala ctic ab sorption by a column density A^ al = 4.35 x 10 2() cm~ 2 
(iKalberla et all : 2005 ). We focused our analysis on the 0.2-10 
keV band of the pn data, and on the 0.3-8 keV band of the MOS 
data. 

We fitted simultaneously the average spectra of the four 
epochs, Dec. 2004, Jun. 2008, Dec. 2008, and 2009 Long Look, 
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0.5 1 2 5 



Observed Frame Energy (keV) 


Fig. 3. Top panel: average 2009 Long Look pn (black), MOS1 
(red), and MOS2 (blue) spectra modelled with a power law emis- 
sion absorbed by the Galactic interstellar medium, along with 
the residuals to the model, scaled by the error bar. Bottom panel: 
spectral residuals to the power law model for the other three 
epochs: Dec. 2008 (green triangles), Jun. 2008 (red circles), and 
Dec. 2004 (blue squares); the data have been rebinned for visual 
purpose only to respectively 5, 5, and 7cr. 


to investigate the causes of the observed variability on month 
time scales. For each epoch, the model parameters were kept 
fixed between the MOS and the pn datasets, except for the pri- 
mary power law normalizations that were left free to vary within 
10% to account for the instrumental cross-calibration uncertain- 
ties. For the sake of clarity, in most of the following plot we show 
the pn data only. 

The spectra were initially fitted to a simple power law model, 
with fixed slope and free amplitude between different epochs. 
Results are shown in Fig. [3] The fit is very poor, with a chi 
square per degrees of freedom ^f 2 /v = 8801/1591, and a very 
flat photon index T ~ 0.6. The spectral residuals of all the four 
epochs have a similar shape, and indeed allowing for the spectral 
slope of the power law to vary between epochs does not improve 
the fit statistics. Three main deviations from the simple power 
law model are consistently observed in spectra from all the three 
instruments and in all the four epochs of observation: (a) both 
negative and positive residuals in the soft band, £<1.5 ke V, (b) 


a deficit of counts at the highest energies probed by the EPIC 
camera, £ ~ 7 - 10 keV, and (c) a cutoff at £ < 2 keV. 

We then tested several models of increasing complexity, in 
order to reproduce the complex spectral shape of PG 1 126-041. 
A neutral absorber fully covering the emission source is not 
required by the data. On the other hand, allowing for the ab- 
sorber to be only partially covering the source (zpcfabs model 
in Xspec) significantly improves the fit statistics with respect to 
the power law model (A^ 2 /Av = 5148/8), but still gives an un- 
acceptable fit (x 2 r ~ 2.3). A better fit for the same number of 
degrees of freedom (A;p 2 /Av = 6567/8) is given by a ionized 
absorber totally covering the source (modeled with XstaiQ), but 
again, the fit is statistically unacceptable Of 2 ~ 1.5). 

The ionized absorber was then allowed to be only partially 
covering the X-ray source (model A), by adding a secondary 
soft power law with the same slope as the primary one. The ra- 
tio between the normalizations of the two power laws gives the 
fraction of the direct (unabsorbed by intrinsic absorption) flux 
£; the absorber covering factor is then C/ = (1 — £). We obtain 
an improvement in the fit statistics of A^/Av = 341/4, and the 
resulting fit is marginally acceptable, with a^ 2 ~ 1.199. 

The model A parameters are reported in Table [3] along with 
their lcr statistical errors. We note how the direct flux fraction 
in every epoch is around 2-3%, as is observed for the scat- 
tered component in absorbed AGN; scattered emission accom- 
panied by recombination emission lines would be naturally as- 
sociated to the presence of ionized gas along the line of sight. 
Unfortunately, with the S/N of our observations we can not dis- 
criminate among the partial covering and scattering scenarios 
in reproducing the soft band residuals, and the secondary soft 
power law physical interpretation is not unambiguous. To the 
first order, model A is able to reproduce the complex spectral 
shape of PG 1126-041, in particular recovering the cutoff at 
£ < 2 keV and most of the soft band spectral residuals. While 
the Jun. 2008 and Dec. 2008 spectra are overall well reproduced, 
the low flux states of the 2009 Long Look and Dec. 2004 still 
show strong residuals at £ < 1 keV and £ > 6 keV. 

In particular, the top panel of Fig. [4] shows a magnified view 
of the 4-10 keV residuals of the 2009 Long Look pn data, rela- 
tive to the simple power law emission (top panel) and to model 
A (lower panel). The two arrows mark strong negative residu- 
als at ~6.7 and 7 keV in the observer frame. We make an initial, 
phenomenological model for these residuals by adding Gaussian 
profiles with negative normalization to model A of the 2009 
Long Look dataset. Even though the model adopted is only a 
rough approximation for a more complex physical context, the 
basic properties such as the FWHM and the EW of the absorp- 
tion features are still physically meaningful. The inclusion of 
two narrow (unresolved, <x s 0 eV) Gaussian absorption lines 
improves the fit statistics by A^ 2 /Av = 60/4 with respect to 
model A. The energies and equivalent widths of the two ab- 
sorption lines are E abs = 7.02+°“ and E abs = 735+0.02 keV, 

EW" bs = -80/ 2 ° and EW“ bs = -127+^j eV, in the source rest- 
frame. Allowing the widths of the absorption lines to be free 
parameters does not improve the fit (A^/Av = 2/2). In the bot- 
tom panel of Fig. [4] we show the 99% confidence contours for 
the absorption lines centroid energy and (negative) normaliza- 


2 We applied the precompiled grid25, publicly available at 

ftp : //legacy . gsfc .nasa . gov/software/plasma_codes/xstar/xspect 
which is computed for an ionizing continuum with T = 2, a gas shell 
with n = 10 12 cnT 3 , a turbulent velocity^,,,,/, = 200 km s _1 , and 
solar abundances, see iKallman & Bautista : d200 1/) and the Xstar 
documentation. 
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Table 3. Spectral analysis results of the simultaneous fit to the 
average spectra of each epoch 



Jun. 2009 

Dec. 2008 

Jun. 2008 

Dec. 2004 


Model A: ionized partial covering 



Xspec : 

zpo 2 + XSTAR*zpoi 


r 


2.15 

+0.01 

-0.01 


N IkeV 

1 05 +015 

i.UJ_oo3 

i oq+0.09 
A -0.08 

1 40 +006 
A -0.06 

1 oo +002 

i.UU_oo2 

Ajm.i. 

14 36+0.61 

ih-.ju_ 047 

4 q i +0.24 
^•^-0.25 

7 Q7+0.36 
' ' -0.24 

9 7Q+0.18 
y ‘ /o -0.08 

log£ mip 

1 66 +002 
i.uu_0 02 

i 44+0.06 
i -^ -0.05 

i co +0.02 
a - jo -0.02 

i C7+O.OI 
A,J/ -0.01 

f'm.i. 

0 97+°- 01 
^• y ' -0.01 

0 98 +0 01 
u -^°-o.oi 

0 98 +0 01 
u -^°-o.oi 

0 97+°- 01 
^• y ' - 0.01 


X 2 /v 1893/1579 


Model B: ionized partial covering + highly ionized absorber 

Xspec : zpo 2 + XSTART v XSTAR2*zpo 1 


r 


2.14 

+0.01 

-0.03 


N IkeV 

1 1 9+O.O6 
1 ■ 1 — 0.13 

1 O9+0.11 

A -° z -o.io 

1 4 1 +0.09 
a, ^ a -0.08 

1 oi +003 

1 ‘^-0.04 

N mi - 

ly W 

14.82$“ 

0 1 7+0.24 
J l -0.12 

6 93+0.57 
u,: * -1.00 

q 7 c +0.24 
y ' ' -0.10 

log 

1-661S! 

1 9Q+0.04 
1 "^ y - 0.04 

1 C4+0.02 
A,J -0.04 

1 C7+0.01 
A,J -0.01 

/-' m . i . 

w 

0.97/™i 

0 98 +0 01 

0 98 +0 01 
u -^°-o.oi 

0 97+ 0 01 
' -0.01 

<'• 

75.18/™= 

< 6.82 

<6.15 

62 54+ 27 - 45 

-35.44 

log 

3 43+0.09 
' , * ho -0.08 

< 3.00 

< 2.73 

9 7C+O.I2 
' -0.20 

V>, out 

0 054 +0 09 

U,UJ -0.03 

0.0 F 

0.0 F 

o.o F 

+ 2 /v 


1771/1570 



Rest Frame Energy (keV) 

4.24 5.3 6.36 7.42 “ 8.48 9.54 10.6 



d, t i l i i i 

'4 5 6 7 8 9 10 


Observed Frame Energy (keV) 



Notes. Units are the following: [10~ 3 ph s~ ! keW 1 cirT 2 ] for the 
power law normalizations Nu^y; [ 1 0 22 cm -2 ] for the column densities 
Nw,N'w and N jj/'; [erg cm s _l ] for the ionization parameters 
and c for the highly ionized absorber outflow velocity v out ; the 
photon index T and the covering fraction C / are adimensional parame- 
ters. Errors quoted are at lcr confidence level, a superscripted F denotes 
a fixed parameter. 


tion, in the source rest-frame. The thin black lines are referred 
to the addition of the absorption lines to the simple power law 
model, while the thick red lines are referred to the addition of 
the absorption lines to model 40 The four vertical green lines 
mark the theoretical energies of the strongest transitions of Fe 
XXV and Fe XXVI, while the two dotted blue lines mark the 
position of the Fe XXV Hen- and Fe XXVI Lycr transitions (rest- 
frame energies E ~ 6.70 and 6.97 keV), both blueshifted by 
Vom = 0.055c. The correspondence with the detected energies is 
striking, so we identify the two absorption lines with these two 
transitions. The moderately ionized absorber (m.i.) responsible 
for the opacity in the soft band has an ionization state too low 
to account for the strong absorption in the iron K band. So we 
added another layer of highly ionized (h.i.) gas modeled again 
with XSTAR. We used the same input parameters as the grid 2 5 
but we find a turbulence velocity v tur b = 1500 km s _I to better 
reproduce the residuals. Only in the case of the 2009 Long Look 
data we allowed the velocity shift of the highly ionized absorber 


3 We note how the EW of the second absorption line is lower when 
added to model A than to the simple power law model. This is because 
in model A some of the high energy negative residuals are accounted 
for by the massive partially covering ionized absorber, in form of neu- 
tral and moderately ionized Fe absorbing photons of < 7 keV in the 
observer frame. 


Fig. 4. Top panel: 4-10 keV spectral residuals of the 2009 Long 
Look pn data, relative to the simple power law model (top panel) 
and to model A (bottom panel). The two arrows mark strong 
negative residuals that are present in both cases. Data have been 
rebinned to 3<x significance. Bottom panel: 99% confidence con- 
tours for the rest-frame normalizations and centroid energy of 
the absorption lines used to model the residuals; the black con- 
tours refer to the addition of the absorption lines to the simple 
power law model, while the red ones to the addition to model 
A. The dotted blue lines refer to the two main spectral transition 
of highly ionized iron, Fe XXV Herr and Fe XXVI Lya, both 
blueshifted by u ou , — 0.055c. 


to be a free parameter; given the much lower S/N of the other 
three epochs, we fixed the velocity shift of this component to 
be null. The addition of the highly ionized absorber significantly 
improves the fit with respect to the moderately ionized partially 
covering absorber model (A^ 2 /Av = 122/9), and overall gives 
a good representation of the data (x 2 r ~ 1.128). Spectral param- 
eters along with their errors are reported in Table [3] model B. 
The corresponding theoretical model is plotted in top panel of 
Fig. [3 while the spectral residuals to the model are plotted in 
the bottom panel of Fig. 0for each epoch. 

Significant residuals at E < 1 keV relative to model B are 
present only in the 2009 Long Look and Dec. 2004 spectra. 
The residuals could be either due to an improper modelization 
of the absorbers, or to an unmodelled component that emerges 
in the soft band when the source is in a low flux state. In any 
case, we know from the timing analysis performed in Sect. [4] 
that in this energy range the source shows no spectral variabil- 
ity over short time scales. We find a statistics improvement of 
A^ 2 /Av = 42/2 and A^f 2 /Av = 23/2 for the addition of two 
narrow Gaussian emission lines with centroid energies and nor- 
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Fig. 5. Top panel: the solid lines show the power law emission 
after the passage through two absorbers, one of which is partially 
covering the source (model B). This model was fit to the aver- 
age 2009 Long Look (black), Dec. 2008 (green), Jun. 2008 (red), 
and Dec. 2004 (blue) epoch spectra. The dashed lines represent 
the direct continuum components that are not absorbed. Bottom 
panel: average 2009 Long Look (black points), Dec. 2008 (green 
triangles), Jun. 2008 (red circles), and Dec. 2004 (blue squares) 
spectral residuals relative to model B; the data have been re- 
binned for visual purpose only to respectively 15, 5, 5, and 10cr. 


malizations fixed between the epochs. The best fit rest-frame en- 
ergies are E\ mi = 0.55+°'° J keV and E e 2 mi = 0.91+°°| keV, while 
the equivalent widths EW™' ~ 10 - 30 eV and EW ™‘ ~ 15-40 
eV, depending on the epoch of observation. The fit statistics after 
the inclusion of the two soft emission lines is^ 2 /v = 1710/1566. 
The observed energies might correspond to O VII K and Fe L 
emission features, but the exact identification is prohibitive with 
the pn spectral resolution. Physically, possible contributors to 
the soft emission can be a scattered component accompanied by 
recombination lines from a photoionized gas, thermal emission 
from an underlying starburst component, or a mixture of the two 
(e.g. IRAS 13197-1627. rMiniutti et al.ll2007h . Unfortunately PG 
1 126-041 is not detected by the RGS, and the pn spectral resolu- 
tion does not allow to distinguish among the different scenarios, 
thus we do not try to further model the soft residuals. 

As for the high-energy residuals, an excess of counts at ob- 
served energies E < 6 keV is present in the 2009 Long Look 


Table 4. 2009 Long Look observation time slices over which 
time resolved analysis was performed. 


Interval 

Exposure 

Count-rate 

/2-10 

/0.2-10 

1 

2 

3 

4 

5 

6 

7 

8 

21.1 

30.0 

5.4 

7.5 

4.6 
5.4 

7.7 
9.2 

0. 180±0.003 
0.156±0.003 
0.212±0.007 
0.195±0.005 
0.215±0.007 
0.194±0.006 
0.338±0.007 
0.245±0.006 

r\ Q'7+0.08 

U ‘ -Q-Q3 

0 71+™ 
' 1 -0.Q3 

1 1 7+0.04 
1 • 1 ‘ -0.05 

0 98 +0 06 

-Q-Q4 

1 1 7+0.04 

098 « 
-Q04 

9 1 1 +0.59 

^Q-o.os 

LOO-- 

0.83/«;| 

1 99+O.O4 
-0.06 
1 19+O.O7 


Notes. Col.(l): Time interval as marked in Fig. [2j Col. (2): Exposure 
time in ks; Col. (3): 0.2-10 keV pn count rate in ct s -1 ; Col. (4): best fit 
pn observed flux computed over the 2-10 keV observed band, in units 
of 10~ 12 erg cnT 2 s -1 ; Col. (5): best fit pn observed flux computed over 
the 0.2-10 keV observed band, in units of 10~ 12 erg cm -2 s~*. Errors 
quoted at lcr confidence level. 


data, suggesting the presence of complexities in the Fe K band. 
We note that at these energies PG 1126-041 is highly variable 
over kilosecond time scales (Sect. ED- The addition of a Fe I 
K a emission line, forced to be narrow and with a centroid en- 
ergy of 6.4 keV in the source rest-frame improves the fit statis- 
tics by A^ 2 /Av = 15/1, corresponding to a significant detection 
(F-test probability > 99.97%). The best fit model statistics is 
X 2 lv= 1691/1565, with a reduced chi square^ 2 = 1.080. There 
are still shallow residuals in the Fe K band, that could be ei- 
ther modeled with a broad Gaussian emission line, a relativistic 
Fe K line, or a partial covering absorber, in all cases suggesting 
the possible presence of further complexities in the iron K band 
(see M. Giustini Ph.D. thesis for details). However, given the low 
statistical significance of the residuals, their appearance only in 
the 2009 Long Look observation, and the high variability of the 
source at those energies over short time scale during the same 
observation, we do not try to model these residuals and turn our 
attention on the time resolved spectral analysis of the 2009 Long 
Look data. 


6. Time resolved spectral analysis 

We perform time resolved spectral analysis on the EPIC pn data 
extracted in the eight time intervals of the 2009 Long Look ob- 
servation marked in the middle panel of Fig. [2] Exposure times 
and background corrected count rates are reported in Table[4]for 
each interval. Given that the count rates are the same in several 
intervals, we merged the 3 rd and the 5 th intervals together and the 
A ,h and the 6 th one together. The unfolded 0.2-10 keV pn spec- 
trum is plotted in Fig. [6] for each of the six resulting intervals. 
At energies £>1.5 keV there are variations of about 50% on 
a few kiloseconds timescales (e.g., from the 2 nd to the 3 rd time 
interval, and from the 1 th to the 8' /! time interval), and variations 
as high as 100% with respect to the average are observed dur- 
ing the 1 th interval. All the spectra were simultaneously fitted to 
the best fit model B described in the previous section, with all 
the parameters tied between the time intervals except than for 
the power law normalizations. The fit to the data is acceptable, 
X 1 /v = 925/889, however the model fails to properly reproduce 
the Fe K band absorption residuals. Allowing for the highly ion- 
ized absorber column density and ionization parameter to vary 
within each temporal slice improves the fit by A^ 2 /Av = 45/6. 
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Energy (keV) 


Fig. 6. The unfolded 0.2-10 keV EPIC-pn spectra extracted in 
the different time intervals marked in Fig. [3 I s ' (black), 2 nd 
(red), 3 rd + 5 th (green), 4 th + 6 th (blue), 1 th (light blue), and 8 th 
(magenta). Data points have been visually rebinned to 10 cr con- 
fidence level. 


Because of the much lower S/N in individual temporal slices, it 
is not possible to assess whether the absorber varied in column 
density, ionization state, or a combination of them. Fluxes com- 
puted using this best-fit model are reported in Table [4] 

We then removed the highly ionized absorber from the 
model, and visually inspected the residuals in each temporal 
slice. Given the similar shape of the residuals, we merged the 
yd + 4 */i + y h + yii t j me intervals together, and the 1 th + 8 th time 
intervals together. The 5-9 keV spectral residuals are shown in 
the left column of Fig.|7]for each of the four resulting intervals; 
time is flowing from top to bottom. There is clear variability 
of the highly ionized absorber over very short time scales: in 
particular, there are no signs of high-energy absorption during 
the 2 nd interval, while the effect of iron absorption during the 
1 th + 8 ,h is dramatic. The right column of Fig.|7]shows the corre- 
sponding 68% and 90% confidence contours for the energy and 
(negative) normalization of a narrow Gaussian absorption line 
used to model the residuals. One can see how the complex ab- 
sorption features at ~ 7 keV in the observer frame do disappear 
within the few ks that elapsed between the 1 st and the 2 nd time 
intervals, and then re-develop during the 3 - 6 th temporal slice. 
The ~ 7 keV absorption trough becomes the deepest and the 
largest during the 7 - 8 lh interval, and its shape becomes quite 
complex, as confirmed by the confidence contours. The intense 
absorption structure at ~ 6.5 keV does develop in a few ks and 
is also visible in the merged 2009 Long Look spectra (top panel 
of Fig. |4}. Summarizing, the continuum power law amplitude 
oscillations dominate the observed spectral variability over ks 
time scales. The highly ionized absorber is also observed to be 
variable, while all the other components (i.e. the photon index, 
the partial covering moderately ionized absorber, and the narrow 
Gaussian emission lines) are found to be constant during the ob- 
servation. 

As for the other epochs, the Jun. and Dec. 2008 observa- 
tions were too short to search for short-term variability. Strong 
spectral variability is found instead during the 2004 observation: 
Fig.[8]shows the unfolded pn spectra extracted during the first 10 
ks (black) and during the last 10 ks (red) of the Dec. 2004 obser- 
vation. The variability pattern is the same as of the 2009 Long 
Look observation i.e. the source is variable on ks time scales 



Observed Frame Energy (keV) Observed Frame Energy (keV) 


Fig. 7. Left panels: 5-9 keV residuals for each temporal slice 
to model B, after the variable highly ionized outflowing ab- 
sorber was removed. From top to bottom: I st (black), 2 nd (red), 
yd + yh + yh + yh (blue), 1 th + 8 th (magenta) interval. Right 
panels: 68% (continuous line) and 90% (dotted line) confidence 
contours for the energy and intensity of a narrow Gaussian ab- 
sorption line used to model the residuals in each slice. 



Energy (keV) 


Fig. 8. The unfolded 0.2-10 keV EPIC-pn spectra extracted dur- 
ing the first 10 ks (black) and during the last 10 ks (red squares) 
of the Dec. 2004 observation. Data points have been visually re- 
binned to 10 cr confidence level. 


only at energies E > 1 keV, and its spectral variabilty is domi- 
nated by continuum amplitude variations. 

7. Optical/X-ray photometry 

Thanks to the presence of the OM aboard XMM -Newton, simul- 
taneous optical/X-ray photometry was performed for each epoch 
of observation. The OM data were all taken in Image mode, and 
were reduced and calibrated using standard pipeline procedures 
through the omi chain SAS task. Background corrected count 
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rates were converted in magnitudes for each exposure (expo- 
sure times are ranging between 1 and 4 ks). We are interested in 
measuring /2500 to compute the optical-to-X-ray spe ctral index 
a m = -0.384 \og(£ 2 50 okl hkev), as defined by Tananbaum et al. 
d!979i) . We used the weighted mean of the magnitudes mea- 
sured in the filters with effective wavelength straddling 2650 A 
in the observed frame (2500 A in the rest frame). The relevant 
filters are the UVW1 and UVM2 (A eff = 2910 and 2310 A, re- 
spectively). We then corrected this magnitude for Galactic ex- 
tinction A2650 = 0.382, estimated using the extinction law of 
ICardelli et al.l (il989i) . Flux densities for each epoch were then 


computed using Vega as a calibrator. Optical photometry results 
are reported in Table 0 

As for the 2 keV flux density measurements, we used our 
best fit model to measure the Galactic absorption corrected rest 
frame fikev for each observation. Flux densities were then con- 
verted to luminosity densities, and the optical-to-X-ray spec- 
tral index was computed. It is well known that the a ox spec- 
tral index in AGN stro ngly depends on the intrinsic UV con- 
tinuum lumin osity ( e.g..mgnali et al 12003b [Strateva et a 1 120051: 
ISteffen et al.l 120061: lYoung et al.l l2010h . We compare the ob- 
served a ox with the value a m (£ 1500 jJ expected on the basis 
of the _2500 A UV luminosity of PG 1126-041, using Eq. 2 
of ISteffen et al.l ( 2006 ). which gives us the parameter Aa^x = 
a ox ~ o; ox ({ 2 5 ooa)- We note that in the UV luminosity range of 
PG 1 126-041 there is a 3<x dispersion in the expected Q' ox (G 5 ooa) 
of about 0.1, so that only values of lAa^xl >0.1 are statistically 
significant (ISteffen et al.ll2006h . We also computed the {^kev ^ u " 
minosity densities corrected for intrinsic X-ray absorption, and 
the corresponding a™ 11 and Ao')’" rr . In Table [6] we report the 
Optical/X-ray photometry results. 

There are variations in the observed luminosity density of 
PG 1126-041 both at UV and at X-ray wavelenghts. While the 
observed 2500 A flux density is steadily increasing from 2004 to 
2009, and is varying by about 30%, the observed 2 keV flux den- 
sity is increasing from 2004 to 2008, and then decreases again 
in 2009, with variations as high as a factor of ~ 4 - 5 between 
the Dec. 2008 and the 2009 Long Look observation. As a re- 
sult, the observed a ox is variable between the different epochs, 
spanning from values typical of radio quiet type 1 AGN (i.e., 
a ox ~ -1.7 in the Dec. 2008 observation) to values typical of 
Soft X-ray Weak AGN (i.e., a ox ~ -2 in the 2009 Long Look 
observation). However, compared to the expected a O x(^2500A)’ 
the source is found to be X-ray weak in all the epochs, with 
Acrox ~ -0.3 - 0.6. Once the effect of intrinsic X-ray absorp- 
tion is taken into account, the maximum observed variations in 
the 2 keV flux density between the different epochs decrease to 
about a factor of two. Consequently, also the Aa^x decreases to 
Aar™ 11 ~ 0.1 - 0.3. It is worth noting that there could be intrinsic 
UV absorption at the redshift of the source, e.g. due by dust in 
the AGN host galaxy; the amount of intrinsic UV absorption is 
unknown so we do not try to model it, but we keep in mind that 
due to this effect the actual (intrinsic) a ox might be flatter than 


8. Discussion 

PG 1 126-041 is observed to be highly variable in X-rays on time 
scales of both months and kiloseconds. The variability on time 
scales of months (observed to be as high as 4x in flux) is domi- 
nated by a spectral component emerging at E < 6 keV (Fig.Q}. 
On the other hand, the RMS variability analysis performed on 
the 2009 Long Look pn observation (92 ks of contiguous good 


exposure time) revealed the presence of a spectral component 
emerging only at E > 1.5 keV that dominates the variability on 
time scales of kiloseconds (up to 50% variations in flux, Fig.0. 

A wealth of information has been obtained both on the intrin- 
sic X-ray continuum emission of PG 1 126-041 and on the repro- 
cessing media that happen to be in the inner regions of this AGN. 
We discuss our results starting from the inner regions around 
the SMBH (i.e., the shortest time scales) and moving farther out 
(i.e., the longest time scales). We adopt the SMBH mass esti- 
mate of Mbh = 1.2 X 10 8 M q _ as given b v iDasvra et al.l d2007l) . 
The Mbh estimated bv IVestergaard & Peterson ( 20061 ) through 
scaling relations is a factor of 2.4 lower. Adopting the bolomet- 
ric luminosity L\ m , = 10 12 L 0 ((Sanders et alll989i) and assuming 
an accretion efficiency 77 = 0.1 gives an accretion rate needed to 
power PG 1126-041 of M ~ 0.7 M 0 /yr and an Eddington ratio 
A = Li,oi / Lb,u ~ 0.26. The gravitational radius is r g ~ 1.8 x 10 13 
cm and the corresponding light travel time is ti ~ 600 s. 

Two absorbers are detected in the X-ray spectra of PG 1 126- 
041, a moderately ionized (log£'"-' ~ 1.5 erg cm s~') and a 
highly ionized (log^ 1 ' ~ 3.5 erg cm s 1 j one. Strong variations 
of both the intrinsic continuum and of the highly ionized out- 
flowing absorber are responsible of the observed kilosecond time 
scale variability. The power law photon index F ~ 2 is found not 
to vary on kilosecond time scales, while the intensity of the in- 
trinsic power law spectrum is following the pattern of the 1.5-10 
keV count rate over the whole observation (see the middle bot- 
tom panel of Fig.[2]and Table|4]i, and doubles in a time interval At 
lasting about 8 ks (the 1 th interval). The variability time scale of 
the intrinsic continuum can set a constrain in the geometrical size 
D of the X-ray emission region through the causality argument, 
D < cAt ~ 2.4 x 10 14 cm ~ 1 3r„. The highly ionized absorber 
is observed to be outflowing at v ou , ~ 16,500 km s' 1 during 
the average 2009 Long Look observation. In the time resolved 
spectral analysis of the same exposure, the outflowing absorber 
is found to be variable over very short time scales (Fig. 0. With 
the much lower S/N of each temporal slice it was not possible 
to assess whether the absober varied in ionization state, column 
density, or blueshift. However, the very short time scale variabil- 
ity suggests that the absorber is very compact and very close to 
the X-ray source, and that we are possibly observing rapid mass 
ejections from the inner regions of the accretion disk, e.g., the 
base of an accretion disk wind. We note in particular how the 
absorption features do develop in a very short time scale, from 
the 6 th to the 1 th interval, together with the strong continuum 
flare. 

By making some assumptions, we can try to estimate the 
mass outflow rate M our associated to the highly ionized absorber. 
In general, the mass outflow rate can be written as: 

M ou , oc A(r)p(r) u out (r ) 


where A(r) is a geometrical factor that accounts on how the 
flow diverges, p(r ) and v out {r) are the density and velocity of 
the flow. By assuming a spherically symmetric, isotropic, steady 
flow with a constant velocity, one can write: 

M ou , = 4 7r r 2 n m H v ou , C f F v ( 1 ) 


where mu is the hydrogen atom mass, r is the absorber distance 
from the central SMBH, n is the absorber number density, C/ = 
cKlj An is the solid angle occupied by the flow as seen by the 
central point source, and Fy is the volume filling factor. 

If we use appropriate values for PG 1 126-041, we can rear- 
range the expression and write 


Mou, = 0.33 (mi) (no) 2 


Uout 

16, 500 km s“ 



F v M q yr 1 
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Table 5. Optical photometry 


Epoch 

V 

[5430 A] 

U 

[3440 A] 

UVW1 
[2910 A] 

UVM2 
[23 10 A] 

UVW2 
[2120 A] 

M(2500) 
[2500 A] 

Dec. 2004 

1 4.7 1 ±0.01 

14.09±0.01 

13.87±0.01 

14.13±0.02 

14.01±0.04 

13.60 

Jun. 2008 

(...) 

(...) 

13.71 ±0.01 

13.92±0.01 

(...) 

13.42 

Dec. 2008 

(...) 

(...) 

13.82±0.01 

13.86±0.02 

(...) 

13.46 

2009 Long Look 

(...) 

(...) 

13.564±0.004 

13.656±0.008 

13.72±0.02 

13.22 


Notes. For each filter we list the effective wavelength; the magnitudes measured in each filter are not corrected for the Galactic extinction, while 
the extrapolated magnitude at 2500 A is corrected for the Galactic extinction. 


Table 6. Optical/X-ray photometry 


Epoch 

(1) 

l°g /2500 
(2) 

log flkeV 

(3) 

log £2500 
(4) 

log £2 keV 

(5) 

(6) 

(7) 

l°g& 

(8) 

<7 

(9) 

A< rr 

GO) 

Dec. 2004 

-25.46 

-30.46 

29.49 

24.49 

-1.92 

-0.52 

25.28 

-1.68 

-0.28 

Jun. 2008 

-25.39 

-30.17 

29.56 

24.77 

-1.84 

-0.43 

25.31 

-1.63 

-0.22 

Dec. 2008 

-25.40 

-29.80 

29.55 

25.15 

-1.69 

-0.28 

25.57 

-1.53 

-0.12 

2009 Long Look 

-25.30 

-30.46 

29.64 

24.49 

-1.98 

-0.56 

25.26 

-1.62 

-0.20 


Notes. Col.(l): Epoch of observation. Col.(2): 2500 A rest frame flux density corrected for Galactic extinction, in erg s 1 cm 2 Hz *. Col. (3): 
2 keV rest frame flux density corrected for Galactic absorption, in erg s _1 cnT 2 Hz -1 . Col. (4): 2500 A rest frame luminosity density corrected 
for Galactic extinction, in erg s _1 Hz -1 . Col. (5): 2 keV rest frame luminosity density corrected for Galactic absorption, in erg s” 1 Hz'. Col. (6): 
Observed optical-to-X-ray spectral index. Col. (7): Difference between the observed a ox and the one expected on the basis of the £2500 luminosity. 
Col.(8): 2 keV rest frame luminosity density corrected for Galactic and intrinsic absorption, in erg s _1 Hz -1 . Col. (9): Optical-to-X-ray spectral 
index corrected for the intrinsic X-ray absorption. Col.(lO): Difference between the corrected Q'“ rr and the one expected on the basis of the £2500 
luminosity. 


where «n = «/ 10 1 1 cnT 3 , rm = r/l0r s . In the expression above, 
the geometrical covering factor 0 < Cf < 1 and volume filling 
factor F v < 1 are unknown, and depend on the wind geometry 
and duty cycle. The density and the radial distance of the flow 
are also unknown. Assuming that the outflowing absorber is a 
thin spherical shell of gas in photoionization equilibrium, one 
can use the definition of the ionization parameter £ = L lo „l nr 2 to 
replace the unknown density and distance with the ratio 
However, in addition to the strong assumptions made regarding 
the spherical geometry of the outflow, there are uncertainties of 
the ionizing continuum luminosity, and especially of the actual 
value of the ionization parameter that may affect our mass out- 
flow rate estimates. In our spectral analysis results we quoted 
only the statistical errors on £: the systematic errors related to 
the input parameters used in the photoionization code (Xstar in 
our case) can be much larger, especially for highly ionized ab- 
sorbers where only a few spectral transitions are available (e.g., 
Fe XXV, Fe XXVI). The major dependence of the best fit £ value 
is related to the slope r, on adopted for the ionizing power law 
continuum. This is because £ oc H r /«H, where n y is the ioniz- 
ing photons number density, ;; H is the photoionized gas number 
density; fixing the ionizing luminosity, a steeper r io „ results in 
a smaller number of high-energy photons able to strip electrons 
from highly ionized Fe than for a flatter r,„„ . This leads to larger 
best fit values measured for £ when using a steeper ionizing con- 
tinuum. To show this dependence, we generated several Xstar 
absorber models using the same input for all the physical pa- 
rameters, but different slopes r,„„ . We then fitted the 2009 Long 
Look average spectra to each of these models, recording the best 
fit £ value measured for the highly ionized absorber responsible 
for the absorption features visible in Fig. [4] Results are plotted 
in Fig. [9] the systematic error on £ are very large, and even re- 



Fig.9. Measured £ for the highly ionized outflowing absorber in 
the 2009 Long Look spectra, with different F , ow used as input in 
the Xstar simulations; Error bars are purely statistical errors at 
lcr confidence level. The vertical dotted line marks the best-fit 
measured F, along with its statistical (dashed lines) and system- 
atic (continuous lines) errors. 


striding the slope of the ionizing continuum to values typical of 
type 1 AGN, 1.8 < F ion < 2.2, gives a systematic uncertainty on 
£ of a factor of ~ 30. The vertical dotted line marks the best-fit 
measured F for the 2009 Long Look spectra, along with its statis- 
tical (dashed lines) and systematic (continuous line) errors that 
have been estimated considering the uncertainties of the power 
law normalization and of the ionized absorbers. The systematic 
uncertainty on is in this case a factor of ~ 6. Just as a com- 
parison with other literature results, we can roughly estimate the 
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kinetic luminosity Lj.,„ associated to the PG 1 126-041 outflow as 
follows. We take the best estimate for g h L from Fig. [9] consid- 
ering the best-fit measured photon index I = 2.15 as the slope 
of the ionizing continuum; for the ionizing luminosity, we take 
the measured unabsorbed 2-10 keV luminosity listed in Table Q] 
Substituting the ratio of these two quantities to the product r 2 n 
in Eq. Q] and using the best-fit outflow velocity, one can write 

Lkin = \M out vl ut ~ 3 (^) x 10 43 erg s- 1 

i.e., less than one hundreth of the bolometric luminosity of PG 
1126-041. 

The uncertainties on the physical quantities related to the 
outflow are huge, and the assumptions made to derive the expres- 
sion for the mass outflow rate and kinetic luminosity are likely 
to be incorrect. The presence of the accretion disk breaks the 
spherical symmetry by itself; furthermore, the highly ionized ab- 
sorber is observed to be variable over very short time scales: the 
geometrical and dynamical effects are most probably very im- 
portant in accretion disk wind scenarios. For these reasons, it is 
very difficult to give a meaningful (i.e. with less than a couple of 
order of magnitudes uncertainty) estimate of the mass outflow 
rate and of its kinetic efficiency. Although one can not exactly 
assess the impact of the fast outflow of PG 1 126- 041 on its en- 
vironment, we stress that its detection is statistically solid and 
model-independent. 

Turning our attention to the larger scales, we note that at 
E < 1.5 keV the power law variability is smeared out by the 
presence of a soft constant component, and at E < 1 keV the 
highly ionized absorber opacity drops (see top panel of Fig. 0). 
so explaining the constant soft flux measured during the whole 
2009 Long Look observation. As for the soft component, we 
modelled it as the direct continuum emission that escapes un- 
absorbed by a partially covering ionized gas. However, the low 
fraction of unabsorbed flux, (1 — C/) ~ 0.03, could physically 
correspond also to a scattered component, that would be natu- 
rally associated to the presence of ionized gas along the line of 
sight. With the moderate pn spectral resolution it is not possible 
to distinguish among the different scenarios, and unfortunately 
PG 1126-041 is not detected by the RGS. However, we stress 
that the adopted modeling for the soft band constant component 
does not influence the absorbers nor the continuum physical pa- 
rameters. 

The PG 1126-041 X-ray variability over long time scales is 
dramatic, and is especially pronounced at energies E < 6 keV. 
The fi kev is observed to be variable up to 4.5 higher in the six 
months elapsed between the Dec. 2008 and the 2009 Long Look 
observations (Fig. Q}. The moderately ionized absorber is de- 
tected in every XMM-Newton observation of PG 1126-041. Its 
column density is significantly varying between the minimum 
N^ 1 ' ~ 3 x 10 22 cm~ 2 measured in the Dec. 2008 observation 
and the maximum N'^; L ~ 1.4 x 10 23 cm -2 measured in the 
2009 Long Look observations. In fact, the moderately ionized 
absorber column density variations dominate the spectral vari- 
ability observed on month time scales, together with amplitude 
variations of the intrinsic continuum up to a factor of two (Table 
0. The ionization state of the moderately ionized absorber 
log <f M - ~ 1.3 - 1.7 is typical of the X-ray warm absorbers, how- 
ever the measured column density are much higher than what 
is usu ally observed in warm absorbers (see e.g. McKerna n et alJ 
l2007h . 

A possible physical scenario i s that of a UV line driven ac- 
cretion disk wind as developed bv lProga & Kallmaiil ( 20041 1 and 


modeled bv ISim et ali (2 ()Tol ). These authors computed the sim- 
ulated X-ray spectra for a non-spherical, hydrodynamical accre- 
tion disk wind in an AGN with Mbh = 10 8 M Q , accreting with 
A = L/LEdd = 0.5, i.e., with conditions not too different to those 
observed in PG 1126-041. The simulated spectra change signif- 
icantly for different inclination angles and different snapshots 
of the flow, confirming both the non-spherical character and 
the highly dynamical behavior of the X-ray spectrum associated 
with such winds. The spectral signature of such flows includes 
both absorption and emission features. For example, the iron K 
band is shaped both by absorption by a highly ionized (Fe XXV, 
Fe XXVI) outflowing phase of the wind, and by emission and 
scattering off the highly ionized base of the flow, that produces a 
Compton bump that mimics a red-skewed iron K line. In the soft 
band (E < 1 keV), there are contributions by a blend of emis- 
sion lines from elements such as C, O, Fe, and bremsstrahlung 
emission by the wind. From a dynamical point of view, the same 
model also predicts a flow component that qualitatively shows 
the same behavior of the highly ionized absorber. At intermedi- 
ate inclination angles above the disk plane and above the “skin” 
of the fast wind, where the gas is most exposed to the X-ray 
ionizing source, recurrent instabilities with blobs (or “puffs”) of 
high density can develop. Some of these puffs quickly become 
overionized and fail to become part of the wind, while others can 
be able to be efficiently accelerated to velocities even higher than 
those of the fast wind, and well above the local escape velocity. 
This flow component would observationally be detectable in the 
X-ray spectra via large column density variations, and sporadic 
high velocity ejections, over short timescales, similar to what ob- 
served during the 2009 Long Look observation of PG 1 126-041 . 
This may suggest that radiation driven accretion disk wind mod- 
els can account in a self-consistent way for most of the observed 
X-ray spectral features in the mini-BAL QSO PG 1126-041. 
Models for the X-ray spectra of accretion disk winds that can 
be fitted to real data are still under developement, therefore we 
must confine ourselves to qualitative considerations. However, 
we note that in this scenario the radial flow assumption is in- 
valid, and the mass outflow rate associated with the highly ion- 
ized phase is overall very different than in the spherical case. 

The moderately ionized absorber, on the other hand, is not 
variable over short time scales but is strongly variable over time 
scales of months. Its observed variability is consistent with the 
scenario suggested to explain the variability of the UV absorp- 
tion lines of the mini-BAL QSO HS 1603+3820 over rest-frame 
time scales of weeks and months dMisawa et al.l2007bll2010l) . In 
this scenario the variability of the UV absorption lines is the re- 
sult of fluctuations of the continuum, which is caused by a screen 
of clumpy, highly-ionized gas between the UV absorber and the 
continuum source. This clumpy screen has similar properties as 
the moderately ionized absorber we have found in PG 1 126-041. 

The observed optical-to-X-ray spectral index is found to be 
highly variable, following the moderately ionized absorber vari- 
ability. Although a ox never gets < -2, when compared with the 
expected a M (^ 25 ooA) based on the UV continuum luminosity it 
is found that PG 1126-041 is observed to be Soft X-ray weak, 
Aa ox = -(0.3 - 0.6) depending on the epoch of observation. 
After correcting for the intrinsic X-ray absorption, the source is 
still slightly X-ray weak, A a c 0 ° rr = -(0.1 - 0.3). The observed 
a ox variability driven by X-ray absorption variability is simi- 
lar to what [Ballo et al. ( 2008 ) observed in the mini-BAL QSO 
PG 1535+547. The 2-10 keV unabsorbed luminosity is still very 
low, (L 2 - 10 ) ~ 2 x 10 43 erg s _1 , and compared with Lboi gives a 
very high bolometric correction, Kboi ~ 200, that however still is 
compatible with the Kboi - a ox relation as in Lusso et al] ( 2010 1. 
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9. Conclusions 

Our XMM-Newton observational campaign on the mini-BAL 
QSO PG 1 126-041 allowed us to characterize the complex high 
energy spectral behavior of the source with unprecedented sen- 
sitivity. We analyzed data for four different pointed observations 
that totally span 4.5 yr. 

The most evident and model-independent result is the detec- 
tion of high column densities of ionized gas along the line of 
sight, that contribute to the observed spectral variability of the 
highly ionized and moderately ionized phase on kilosecond and 
month time scales, respectively. In particular, the highly ionized 
absorber is found to be variable on very short (a few ks) time 
scales, suggesting an origin very close to the SMBH. Also the 
intrinsic X-ray continuum emission is observed to be variable in 
intensity, over both short and long time scales. 

Overall our findings are qualitatively consistent with 
radiation-driven accretion disk wind models predictions, where 
one expects a high column density of X-ray absorbing gas 
shielding the portion of the wind that is accelerated by UV pho- 
tons. The hot X-ray component of the wind is also expected to 
be highly variable with time in such a model. The mini-BAL 
appearance of the wind of PG 1126-041 is qualitatively consis- 
tent with an orientation effect, in which our line of sight and the 
plane of the accretion disk make a larger angle than when view- 
ing classical BAL QSOs. We are looking at a narrower range of 
UV velocities through the wind (hence the small width of the 
UV absorption features) and at the transition zone between the 
fast UV wind and the hot X-ray component of the flow, where 
we expect strong X-ray spectral variability. 

The present X-ray observational campaign of the mini-BAL 
QSO PG 1126-041 has demonstrated that long and short term 
variability studies are quite powerful tools for identifying the 
physical mechanisms at work, and for mapping the dynamics of 
the inner accretion/ejection flow in AGN. 
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